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ABSTRACT 

We present the R-band spectral line lists adopted by the Apache Point Observatory 
Galactic Evolution Experiment (APOGEE). The APOGEE line lists comprise astro- 
physical, theoretical, and laboratory sources from the literature, as well as newly eval¬ 
uated astrophysical oscillator strengths and damping parameters. We discuss the con¬ 
struction of the APOGEE line list, which is one of the critical inputs for the APOGEE 
Stellar Parameters and Chemical Abundances Pipeline, and present three different ver¬ 
sions that have been used at various stages of the project. The methodology for the 
newly calculated astrophysical line lists is reviewed. The largest of these three line lists 
contains 134,457 molecular and atomic transitions. In addition to the format adopted 
to store the data, the line lists are available in MOOG, Synspec and Turbospectrum 
formats. The limitations of the line lists along with guidance for its use on different 
spectral types are discussed. We also present a list of R-band spectral features that are 
either poorly represented or completely missing in our line list. This list is based on the 
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average of a large number of spectral fit residuals for APOGEE observations spanning 
a wide range of stellar parameters. 

Subject headings: keywords: Astronomical Instrumentation, Methods and Techniques: 
methods: laboratory: atomic - Physical Data and Processes: atomic data - Physical 
Data and Processes: line: identification - Physical Data and Processes: molecular data 


1. Introduction 


The Apache Point Observatory Galactic Evolution Experiment (A POGEE) is one of the pro¬ 
grams that was carried out on the Sloan Foundation 2.5-m Telescope (I Gunn et al.ll2006l ) by the 
third stage of the Sloan Digital Sky Survey (SDSS-III; Eisenstein et al. 2011 ). APOGEE obtained 
high resolution (R ~ 22,500) and high signal-to-noise rat io (S/N > 100) spec tra in the i^-band 


(1.51-1.70 fim) for more than 100,000 cool giant stars (see IZasowski et al 


2013, for more informa¬ 


tion about targeting) spanning all components of the Milky Way (jMaiewski et al.ll2015l ). Stellar 
parameters and indi vidual chemical abundances are derived from the combined APOGEE spectra 
(jNidever et al.ll2015l ) with the APOGEE St ellar Parameters and Che mical Abundances Pipeline 
(ASPCAP), which is described in detail in Garcia Perez et al. ( 2015 ). ASPCAP uses a grid of 
synthetic spectra to determine stellar parameters and abundances by finding the best match to the 
observed spectra interpolating within the grid. For the public da ta releases, the s ynthetic spectra 
have been calculated using the spectral synt hesis code ASSeT (IKoesterkel l2009l ). which itself is 
based, in part, on the synthesis code Synspec (jHubenv fc Lanal201lh . 


In order to run ASPCAP on the first year of APOGEE results (DRIO; lAhn et al.ll2014l) . an 
initial line list was generated and adopted. In subsequent years, improvements were made and 
different methodologies were adopted in th e line list used for the release of the full three year data 
set of APOGEE (DR12; Alam et al. 2015 ). The APOGEE internal naming scheme for the DRIO 
and DR12 line lists are 20 1 10510 and 2(313121 6, res pectively. In t he ye ars between DRIO and DR12, 
Garcia Perez et al.l (|20I3l ). ISmith et al.l (l2013l ) and ICunha et al.l (j2015l ) made use of an intermediate 
line list, 20120216. The APOGEE naming scheme includes the year, month and day label to keep 
track of changes made to each line list. Because these dated names are long and difficult to associate 
with each data product we will refer hereafter to 20110510, 20120216 and 20131216 as DRIO, INT 
and DR12 line lists, respectively. The line list continues to evolve in APOGEE-2, an extension of 
the project as part of SDSS-IV. 


In this paper we document the methodologies employed by APOGEE to produce the three 
H-hand line lists that were used by ASPCAP (and several independent analyses) to derive stellar 
parameters and chemical abundances. The first section of this paper describes the base line list 
taken from literature sources, including those derived from laboratory, theoretical, and astrophysical 
sources. Parameters of the transitions, even if they are well studied, have associated uncertainties. 
We aspire to improve on the theoretical line parameters by comparing synthetic spectra with 
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the observed high-resolution spectra of two well known stars (the Sun and Arcturus) to produce 
“astrophysical” oscillator strength values and damping constants. In the second section, we describe 
the code used to derive our astrophysical log{gf) and damping parameters. We also describe some 
issues identified in the INT and DR12 APOGEE line lists and describe the impact these issues 
have on the stellar parameters and abundances derived from those lists. In Section 0] we detail 
the stellar features which appear to be missing from our line list, based on synthesis from a large 
portion of the DR12 stellar library. In the last section, we describe the line list formats and review 
the performance of the line list as described in a number of papers in the literature. 


2. Base Line List from Literature 
2.1. Oscillator Strengths 

This section describes the various literature sources that were considered for the base line list. 
Criteria for accepting and rejecting various sources are presented. When we refer to the oscillator 
strengths, we use the standard expression log(g'/) for the base ten logarithm of the product of the 
lower level degeneracy and absorption oscillator strength. 


2.1.1. Molecules 

The molecular line lists are taken almost entirely from literature sources. Eor DRIO an attempt 
was made to fit the molecular features to the very weak lines seen in the Sun in order to define 
solar gf values. Eor the INT and DR12 line lists, the molecular astrophysical log(g'/) changes were 
removed in favor of adopting the best literature values available. Below we discuss separately each 
molecule and note references adopted, as well as other line lists that were tested. 


2.1.1.1. CN We used the Kurucz (CD-ROM 18) CNAX.A SC and CNBX.ASC l i sts as the base 
line list. This combined list was tested against the line list in iMelendez &: Barbuvi (|l999l . MB99). 
The latter is a hybrid of theoretical log(qf ) values provid e d by S. P. Davis and astrophysical 
values fit to the Sun using abundances from ICrevesse et al.l (jl996|) . Our tests indicated that the 
line positions and strengths from MB99 were found to overall provide a better match to the CN 
features in the Sun and Arcturus. Thus, the MB99 list was adopted, while we kept in our line list 
any of the CN lines which were in the Kurucz line list and not in MB99. We tested this composite 
list against the 2010 version of the Plez CN line list (private communication) and found that 
our composite list better fit the Sun and Arcturus spectra, i.e., resulting in fewer poorly fit lines 
and fewer predicted lines without associated observable features. Eor DRIO we changed the CN 
molecules’ \og{gf) values by -1-0.03 to fit the strongest lines in the Sun based on a visual inspection. 
Eor INT and DR12 we removed the astrophysical \og{gf) offset and the composite line list (MB99 
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+ Kurucz) was adopted in the final line list. 


2.1.1.2. CO We tested the Kurucz (CD-R OM 18) COAX.ASC and COXX.ASC composite line 
list against the line list from iGoorvitchl (jl994l l and found the latter to be slightly better in strength, 
but not in wavelength. Thus we adopt the Goorvitch \og{gf) values for the CO lines in the Kurucz 
list, while keeping the original wavelengths as in the Kurucz line list.. There were no lines in the 
Goorvitch list that were not in the Kurucz list. We retained all lines in the Kurucz list that were 
not included in the Goorvitch list. For DRIO we changed the CO molecules \og{gf) values by -0.10 
to fit the strongest lines in the Sun based on an overall visual inspection of the obtained fits. For 
INT and DR12 we removed any astrophysical \og{gf) values and the Goorvitch and Kurucz hybrid 
line list was adopted. 


2.1.1 .3. OH We tested the Kurucz (CD-ROM 18) OH.ASC list against that of I Goldman et al 


(|l998ll and found t he Goldman line s to pr ovide a better ht to the Sun and Arcturus. We adopted 
the OH lines from iGoldman et al.l (jl998l l for our base line list. For DRIO we changed the OH 
molecules log{gf) values by -0.07 to fit the strongest lines in the Sun based on an overall visual 
inspection of the obtained hts. For INT and DR12 we removed any astrophysical log( 5 '/) values, 
and the Goldman-|-Kurucz OH lines were adopted as the final product. 


2.1.1. 4. H 2 The Kurucz (CD-ROM 18) H2.ASC was adopted as a base line list but was sup¬ 
plemented by a few additional lines from hdxx.asc (Kurucz web site, 

http://kurucz.harvard.edu/molecules/). For DR12 we reduced the \og{gf) value for the H 2 
line at 16586.664 A by 6.0 dex, based on a poor fit to the Arcturus spectrum and to be in line with 
the other log{gf) values for the H 2 lines. 


2.1.1.5. C 2 The Kurucz (CD-ROM 18) C2AX.ASC, C2BA.ASC, C2DA.ASC, and C2EA.ASC 
files were used as a base line list. For DR12 we hxed some C 2 AX fe atures that were cl early 
discrepant in Arctu rus by replacing their log(g'/) values with those from iBrault et al.l (jl982l l and 


Kokkin et al.l ( 20071 1. 


2.1.1.6. SiH We used the Kurucz (CD-ROM 18) HYDRIDES.ASC list for the SiH features as 
a base line list. While these features should be visible in only very few giant and dwarf stars, they 
have been included in the line list for completeness. 


2.1.1 .7. FeH For DR12, several FeH line lists were tested, including lines fromlLanghoff &: Bauschlicher 
(|l990ll as implemented in the Uppsala spectral synthesis code BSYN, and iDulick et al.l (|2003ll as 
implemented by Kurucz in fehfx.asc (Kurucz web site, 

http://kurucz.harvard.edu/molecules/). We tested these on H-band NIRSPEC spectra of the 
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cool dwarf stars GL436 (jPrato et alJl2002l l and GL763 (jBender et alJl2005ll . kindly provided by G. 
Bender. Based on the synthetic spectra generated using these two line lists and stellar parameters, 
we expected to find a large number of weak FeH features. We searched for conspicuous spectral 
features in the observed spectra at locations that both line lists predicted there to be relatively 
isolated FeH lines, but did not find any consistently accurate predictions. In addition, we tried to 
determine if FeH lines would help in a statistical sense even if the strongest features could not be 
identified. The hrst test was a cross-correlation of a synthetic spectrum of pure FeH against the 
observed spectra; no peak was found at the rest velocity of either star. The second test was to 
subtract the synthetic spectra with and without the FeH features from the observed spectrum; in 
both cases the scatter increased with the inclusion of the FeH lines. From these tests we conclude 
that the current FeH line lists would not assist in our analysis of the APOGEE spectra, thus the 
EeH lines were not included for DR12. 


2.1.2. Atoms 

Eor the atomic features, we compiled an atomic line list from a variety of different sources. As 
a base line list, we started with the Kurucz line list gfhy3000.dat (Kurucz web site, 
http://kurucz.harvard.edu/linelists/gfhyperlOO/). In this line list we included, in sepa¬ 
rate columns, lab data, the “best” empirical log(g'/) values in the literature, and our astrophysical 
log(fl'/) values. This line list contains many more lines than are typically detectable in //-band stel¬ 
lar spectra of cool giants, but all lines were retained in the line list as it may aid future investigations 
of extremely hot stars, investigations of nebular features, or laboratory efforts. 


2.1.2.1. Laboratory Data The gold standard for line list data is high quality laboratory 
measurements of energy levels, branching fractions and lifetimes. With this type of data we will 
have not only the best quality values but also quantifiable errors that can be used to constrain 
empirical changes from laboratory values, as will be discussed in Section O In this section we 
consider both true laboratory measurements and any theoretical measurements that have such well 
constrained theoretical uncertainty as to allow uncertainties to be included. For example, in the 
case of hydrogen (Paschen lines), the National Institute of Standards and Technology database 
fNIST: iKramida et al.ll2014l ) gives the uncertainties in the log(g'/) as AAA or < 0.3%. Within this 
paper we will refer to any log( 5 f/) that has quantifiable errors as laboratory data. Tabled] contains 
the origins of the different lab sources adopted in our list. Most of these were compiled in the NIST 
database and we adopt their us ual grade-to-uncert ainty conversion (Table ID). In the case o f Ti I, 
we adopted log( 5 f/) values from iLawler et al.l (|2013l l over Blackwell-Whitehead et al.l (j200m when 
available. 


In addition to the laboratory log{gf) values from the sources above, for the DR12 line list we 
updated the theoretical wavelengths with wavelengths from the following sources: 
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Table 1: Laboratory Sources for Oscillator Strengths 


Species 

Source 

H I 

NIST - Wiese V Fuhr 120091 

He I 

NIST - Wiese V Fihr 120091 

C I 

NIST ~ Wiese et al. 119961. NIST - Wiese V Fuhr 120071 

C II 

NIST - Wiese et al. 119961. NIST - Wiese V Fuhr 120071 

C III 

NIST ~ Wiese et al. 119961 

C IV 

NIST - Wiese et al. 119961 

N I 

NIST - Wiese et al. 119961. NIST - Wiese V Fuhr 120071 

N II 

NIST - Wiese et al. 119961. NIST - Wiese V Fuhr 120071 

N III 

NIST - Wiese et al. 119961 

N V 

NIST - Wiese et al. 119961 

0 I 

NIST ~ Wiese et al. 119961 

0 II 

NIST - Wiese et al. 119961 

0 III 

NIST - Wiese et al. 119961 

Na I 

NIST - Kellerher &: Podobedova 120081. NIST - Sansonetti 120081 

Mg I 

NIST - Kellerher & Podobedova 120081 

Mg II 

NIST - Kellerher & Podobedova 120081 

A1 I 

NIST - Kellerher & Podobedova 120081 

Si I 

NIST - Kellerher & Podobedova 120081 

Si II 

NIST - Kellerher & Podobedova 120081 

Si III 

NIST - Kellerher V Podobedova 120081 

S I 

NIST - Podobedova. Kelleher k, Wiese 120091 

Ar I 

NIST - Wiese. Smith k Miles 119691 

K I 

NIST - Sansonetti 120081. NIST - Wiese. Smith k Miles 119691 

Ti I 

Lawler et al. 120131. NIST - Blackwell-Whitehead et al. 120061 

Ti II 

Wood et al. 120131 

V I 

NIST - Sansonetti 120081. NIST - Wiese. Smith k Miles 119691 

Fe I 

Ruffoni et al. 120131 


The species with NIST indicated as a Source contain the noted references in the NIST database. 
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Ti I wavelengths from ISalomanI (|2012l l 


Ti II wavelengths from ISalomanI (|2012l l 


V I wavelengths from iThorne et alJ (|201l|) 


Cr I wavelengths from ISalomanI (j2012l l 


Cr II wavelengths from ISalomanI (j2012l l 


Rb I wavelengths from ISansonettil (|2006l l 


2.1.2.2. Literature Astrophysical and Theoretical log{gf) Values The goal of these ef¬ 
forts is for the line list to be a comprehensive list of all 77-band transitions which may appear in 
APOGEE spectra. With this goal in mind, we have taken the base line list (as defined in 12.1.21) 
and augmented it with additional lines from various theoretical predictions. Since we are not the 
first to generate i7-band line lists, we also tested some of these literature compilations against our 
base list to determine if we might improve the list with transitions from these sources. Since there 
are few 77-band transitions known for elements heavier than copper, we placed extra emphasis and 
effort on Z > 30 transitions. 


A careful review of all pages of the Kurucz web site (http: //kurucz. harvard. edu/atoms. html) 
revealed a few ad ditional Ti, F e , and Cr tr ansitions, w h ich w ere added to the base line list. We 


tested the MB99, iRvde et al. (1200911 . and iRvde et al.l ( 2010l i line lists and found their log{gf) 


values to provide a better fit to the Sun than the base line list. This should not be surprising as 
these three references are based on astrophysical solar log( 5 '/) values. All of thes e log{gf) values 
were added, giving preference to the most recent line lists, e.g.. lRvde et al. (1200911 over MB99 , for 
any features in common between the lists. We added Ce III line s from lWvart &: Palmer il (jl998l ) and 
Y II lines from the DREAM compilation (iBiemont et al.lll999l . Feb. 2011 download ). We replaced 
all log (q f) values for Ca I and Ca II with the values given in iHansen et al.l (Il999l b and ILaughlin 
(j 19921 ). respectively. 


2.I.2.3. Hy p >er-f ine Splitting For DR12 we adopted the energy levels for V I lines from 


Thorne et al 


T 

M 


201 Ih. and we added the V hyper-f ine sp litting (HFS) componen ts using the coeffi¬ 


cients of lPalmeri et al.l (ll997lhlPalmeri et al.l (jl995l i. and lGiizelcimen et al.l (|201llb iGiizelcimen et al 
(j2014l ) and IWood et al.l (120141 ) have published laboratory measurements for the V energy levels since 
the DR12 line list was generated. These values produce HFS components that are i n excellent agree¬ 


ment w ith the earlier theoretical predictions. Based on the Na I components from (j Safronova et al 


(j 19991 ). iHapperl (jl974l )l the HFS is too small to impact the line profile or abundance analysis. 


For Al the HFS can be significant an d its exclusion from the APOGEE line li sts deserves com - 
ment. Using the HFS components from Falkenburg fc ZimmermannI (|l979l l and ISur et al.l (j2005l i 
we find that for the weaker Al lines the impact of HFS is less than 0.1 dex (a typical uncertainty 














































































































for APOGEE) while for the strongest line at 16755 it has offsets as large as 0.4 dex for the coolest, 
most metal-rich giants. ASPCAP measures abundances from all of the A1 features and thus will po¬ 
tentially have a bias toward abundances which are slightly too high and a function of line strength 
(largely metallicity and effective temperature). Examination of the slope and s catter of the A1 


abund ances for open clusters in the APOGEE calibration sample (Eigure 5 in iHoltzman et al 


(j2ni5l ) suggests that this bias is less than 0.1 dex. A secondary check can be made by compar¬ 
ing litera ture determinat i ons of Al in globular clusters compared to those using an APOGEE line 
list as in iMeszaros et al.l (|2015l b Eigure 4. If HES was a significant influence in these stars then 
one would expect that the [Al/Ee] enhanced abundances would be biased positively compared to 
the [Al/Ee] ’’normal” stars and that the affect would be strongest in the most metal-rich clusters. 
There is a small trend among the clusters of the order of 0.2 dex. While Al HES will be added in 
future versions of the APOGEE line list, users of this line list in the future sh ould be aware of th is 
limitation and perhaps avoid the use of the strongest Al feature as was done in ISmith et al.l (120131 1. 


2.2. Damping 


The line width for a spectral feature is a complicated function of stellar rotation, thermal 
broadening, turbulence and surface convection, and several other types of quantum mechanical 
broadening that mostly impact strong lines, including Stark, resonance, and van der Waals. Broad¬ 
ening coefficients are often included in line lists, and one of the most important of these is the van 
der Waals broadening. In this work, what we refer to as “damping” is actually the log of the broad¬ 
ening coefficient: van der Waals collisional damping divided by the number density of hydrogen, or 
logirQ/Nu) = where V is the velocity (set by thermal motions with Tg// = 10000) and 

Cq is the interaction constant. There are only a few different sourc es for the damping or C r values 
in the i7-band: the Kurucz compilation, MB99, Rvde et al. (2009), and Rvde et al. ( 2010 ). MB99 
ac tually used several different r netho ds to get the damping: fits to the solar spectrum, values out 
of iBarklem. Anstee fc O’Maral (jl998l ) and other references. 


The Barklem web site (http: //www. astro .uu. se/$\sim$barkleni/, v2.0 iBarklem. Anstee &: O’Mara 


1998|) has codes for calculating the van der Waals damping. However, these codes only work for 
certain values of the effective principal quantum number. We downloaded v2 and ran it for all 
of the transitions within range. Eor many of the IR transitions the effective principal quantum 
number is outside of the Barklem grid. 

It is possible to approximate the Cq values as 6.46 E—34Ar^, wh ere Ar is the unit-less difference 
in the mean square radius of the two energy levels (|Unsoldlll955h . The radius can be described 
roughly as r = 0.5 * * (5 * -|- 1 — 3 * ((( -|- 1)) (MB99) where n = 1/{^EP/x)^'^, AEP is 

the difference between the energy level of the transition and the ionization energy (y), and I is the 
orbital angular moment um quantum n umber. This approximation yields a Cq value 10 times larger 
than that calculated by lUnsoldl (1195, 5l ) . 
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We assumed the Barklem values to be the best possible source for theoretically derived damping 
constants. In Section [3^ we show that these theoretical values are very close to the calculated values 
based on comparisons to the solar spectrum. Figure [T] shows the difference between the different 
damping constants and the Barklem values as a function of atomic mass number for those lines in 
common. The Kurucz damping constants are systematically too small by several tenths of a dex, 
while the Cq approximation and MB99 values are too large by a few tenths of a dex. From this 
assessment, our ranked preferences for adopted damping values for the base line list are: 


1 . 

2 . 

3. 

4. 

5. 


Barklem. Anstee fc O’Maral (jl998l ) whenever possible. 


Rvde et al. (2009|) and iRvde et al.l ( 20ld i. 


MB99. 

The Cg approximation given in this section. 
Kurucz values 


3. Semi-automated “Astrophysical” Line List Calculations 

It is possible to improve upon theoretical log{gf) and damping values by comparing line profiles 
and strengths of well understood spectra to synthetic spectra. Even laboratory measurements with 
significant error bars could be improved by this type of comparison if done within the errors of 
the measurement. We refer to these corrections based on observed spectra as “astrophysical” line 
parameters. In this section we describe the software and methodology we created to generate 
astrophysical line parameters and compare those parameters to the base line list values. It should 
be noted that in making changes to the atomic line parameters to match the observations one is 
also masking systematic errors in the model atmospheres and the line formation calculations. 


3.1. Astrophysical Software and Methodology 


We have developed a code that can vary the log{gf) and damping values in a line list to match 
one or more observed spectra. The final product of the code is a set of astrophysical log{gf) and 
damping values that we adopt in the final line list. This code relies on having accurate stellar 
parameters for the observed stellar spectra and a constrained and accurate base line list. The code 
was developed within the APOGEE team by DB and is available through the Astrophysics Source 
Code Library (iBizvaev fc Shetronell2ni5l ). This code has evolved, and the version used for DR12 is 
more complex than that used for DRIO. Below we document the state of the code as implemented 
for DR12 and describe the significant changes that were made after DRIO. 
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Fig. 1.— Difference between the log damping values from various sources and the Barklem values. 
The X axis shows atomic number. The red squares are MB99 - Barklem, the blue triangles are Ug 
approximation - Barklem, and the green circles are Kurucz - Barklem. The point type has been 
offset slightly to avoid them overlapping each other. The dashed line at A = 0 is provided to guide 
the eye. 
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The code is written in IDL as a wrapper around the spectral synthesis code MOOG (jSneden 


1974). For DRIP the 2009 version of MOOG was used, and for DR12 the 2013 version was used. 


Zamora et al.l (j2015l i conducted an analysis of the differences between synthetic spectra computed 
with both MOOG and ASSeT and found the differences to be very minor for the Sun and Arcturus, 
< 0.01 except around the hydrogen lines. The astrophysical ^'/-fitting code starts with an input 
line list, which in our case is the best laboratory and literature values available between 1.5 and 1.7 
/im (see Section HL For spectral com parison we adopt an i7-bar id center-of-disk ( COD) atlas for 
the Sun (jLivingston &: Wallace! 1 199 il l and an Arcturus flux atlas (jHinkle et al.lll995l b The adopted 
solar and Arcturus stellar parameters and abundances are listed in Tables [S] and [H The DRIO solar 
model was created from a grid of Kurucz Atlas models (jGastelli et al.lll997ll with computed Opacity 
Distribution Functions (PDFs). We interpolate within this grid using the code “makekuruczpublic” 
(jMcWilliam et al.lll995l b Unfortunately, after the DRIO calculations were made, a mistake entered 
into our methodology and the synthesis of the Sun was done as if it were a full flux spectrum, 
rather than a GOD calculation. This is noted in Table [3] as “not GOD”. In Section [3. 2. II we discuss 
the impact of this on the INT and DR12 line lists and how it may impact the subsequently derived 
abundances. 


2003 


The DR IO Arcturus model atmosphere was interpolated from the 2005 MARGS grid (jGustafsson et al, 
20081 . , further expanded by B. Edvardsson 2015, private communication). The MARGS 


models have opacity sampling instead of ODFs; however, the MARGS models are spherical, with 
appropriate a-enhancement at the metallicity of Arcturus, and thus were deemed to be more ap¬ 
propriate for DRIO. Although the differences between the most re cent Kurucz and MA RGS models 
are minimal at the effective temperature and gravity of Arcturus (IZamora et al.ll2015l ). we adopted 
a Kurucz model for Arcturus after DRIO to be consistent with the rest of the ASPGAP models. 
For the INT line list, we use the same Kurucz grid from which we pulled the solar model. We 
should note that the exact details of the Arcturus model atmosphere for DRIO and INT will have 
little impact on the final line list because, as will be described in more detail below, the final solar 
log( 5 f/) fitting removed all but the weakest or most temperature/gravity sensitive features in the 
Arcturus log(g'/) fitting. For DR12 the log(g'/) htting methodology changed so the model was far 
more critical. For DR12 we adopted a model generated in the same way as the mo del atmosphere 


grid p oints were generated for ASPGAP, namely from the Kurucz Atlas 9 code Isee iMeszaros et al 

201 2I 1. 


The solar abundances for DRIO an d INT were adopted from lAsolund et al.l (|2009l l. while for 
DR12 we adopted the abundances from lAsolund et al.l (|2005l l. This change was made to make the 
assumed abundances in t he line list consistent with those adopted in the Kurucz model atmospheres 
calculated for APOGEE ( Zamora et al. 2015 : Meszaros et al. 2012). 


The Arcturus abundances were set with a number of considerations. Since we were concerned 
about mishtting t he atomic features ii ear molecular features, we forced the G, N, and O abundances 
to match those of ISmith et al.l (120131 ) . who derived those abundance s in an independen t analysis 
in the 77-band using line list INT. The G, N, and O abundances from ISmith et al.l (j2013l l are very 


























































Table 2; Adopted Uncertainties for the Oscillator Strengths 


Grade 

Uncertainty 

AAA 

< 0.3% 

AA 

< 1% 

A+ 

<2% 

A 

< 3% 

B+ 

< 7% 

B 

< 10% 

c+ 

< 18% 

c 

< 25% 

D+ 

< 40% 

D 

< 50% 

E 

> 50%. 


Grade-to-uncertainty conversion for the oscillator strengths as defined by NIST 


Table 3; Model Atmospheres Adopted for Astrophysical log{gf) Value Calculations 


Model 


log(g) 

[Fe/H] 

Vt 

Notes 

Sun 

BRIO 

5780 

4.40 

+0.00 

1.10 

Kurucz OBFNEW 

INT 

5777 

4.44 

+0.00 

1.10 

Kurucz OBFNEW; not COB 

DR12 

5777 

4.44 

+0.00 

1.10 

Kurucz OBFNEW; not COB 

Arcturus 

BRIO 

4290 



1.67 

MARCS 

INT 

4286 

■n 

■g 

1.70 

Kurucz NOVER 

BR12 

4286 

■n 

-0.52 

1.70 

new revised Kurucz 


Note. — Not COD: these syntheses were mistakenly calculated as full flux calculations within MOOG and not set 
as center of disk (COD). 
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similar (within ~0.05 - 0.10 de x) to those derived r ecently bv lAbia et al.l (j2012l l. who used NIR lines 
of CO, OH, and CN, as well as ISneden et al.l (j2014l l. who used optical C 2 , [O I], and CN transitions. 
The mean and standard deviations of these 3 studies are A(C)=8.01±0.05, A(N)=7.66±0.02, and 
A(0)=8.64±0.03, with ^^C/^^C=7.4±1.4; the Arcturus C, N, and O abundances are therefore 
well-constrained in these recent independent studies. We wanted the abundances to be largely 
self-consistent with those adopted in the model atmosphere, so the abundances of Mg, Si, S, and 
Ti were all given a-enhanced values of -1-0.4 above scaled solar, whereas all of the odd Z elements 
(except Al) and iron peak elements were given scaled solar values. A1 has several very strong 
features with several other important features in it s wings, thus we ad opted a value to fit all three 
strong lines. This adopted value was close to that of ISmith et al.l (j2013l l. who excluded the strongest 
line from their analysis. The abundance of C a clearly deviates s ignificantly from the other alpha 
elements, so we adopted a value close to that of ISmith et al.l (j2013ll . For DRIO the Arcturus log{gf) 
values were only adopted for a small subset of those lines-of-interest (LOI) which were above the 
threshold in Arcturus but not in the Sun. For DR12 the adopted Arcturus abundances play a far 
more significant role because the final Arcturus and solar log{gf) values were averaged together, 
although the solar values were given more weight. In retrospect, the self-consistency requirement, 
i.e., forcing the abundances to match the model, should not have been a significant driver and may 
have resulted in a less optimal l ine list that resulted i n an abundance inconsistency for Arcturus in 
ASPCAP (see Section 5.5.1 in Holtzman et al. 2015 1. 


In order to match the line broadening in the Arcturus atlas spectrum, we convolved all synthetic 


spectra with an inst rumental profile and a rotational profile (usini = 2.0 km s (jGray 


1981 


Gray fc BrownI 120061 ). The limb darkening for Arcturus in the NIR is assumed to be 0.46 (iGlaret 
2 OOOI ). Solar synthetic spectra were compared with the COD atlas spectrum, and therefore were not 
corrected for rotation but they were convolved with an instrumental profile. The instrument profile 
was determined from an analysis using cross-correlation against a synthetic spectrum generated 
from the base line list but avoiding the strongest lines. 


For DRIO and INT we selected the LOI using the line strength (opacity at line center) tabu¬ 
lation built into MOOG which allowed us to remove extremely weak lines. For DR12 we refined 
this by running each line individually and noting the line depth. The number of strong lines in a 
spectrum varies significantly with stellar spectral type and with an adopted cutoff for the minimum 
line strength. We calculated the maximum line depth with respect to the normalized continuum for 
each line in the initial line list. Table [5] shows the number of atomic and molecular features deeper 
than the listed depth in the solar and Arcturus spectra. Since the target S/N ratio for APOGEE 
program stars is 100, we considered evaluating only those lines that were deeper than 0.001 with 
respect to the normalized continuum. For reference the weakest lines visible in Figure [3] have depth 
of 0.004 and equivalent widths of <0.6mA and are well below the detection threshold in APOGEE 
spectra. The final set of LOI is the joint list for the Sun and Arcturus corresponding to the >0.001 
level (see Table [5]). Future improvements to the line fitting code will likely involve the entire line 
strength and not the depth at the line center as this would remove line depth variations due to 


































Table 4; Abundances Adopted for Astrophysical log{gf) Value Calculation 


Atomic Sun Arcturus 


Number 

Species 

DRIO 

INT 

DR12 

BRIO 

INT 

DR12 

6 

C 

8.43 

8.43 

8.39 

7.88 

8.27 

7.96 

7 

N 

7.83 

7.83 

7.78 

7.21 

7.26 

7.64 

8 

0 

8.69 

8.69 

8.66 

7.59 

8.54 

8.64 

11 

Na 

6.27 

6.27 

6.17 

5.73 

5.65 

5.65 

12 

Mg 

7.53 

7.53 

7.53 

7.24 

7.41 

7.41 

13 

A1 

6.43 

6.43 

6.37 

6.24 

6.15 

6.15 

14 

Si 

7.51 

7.51 

7.51 

7.02 

7.39 

7.39 

16 

S 

7.15 

7.15 

7.14 

6.70 

7.02 

7.02 

19 

K 

5.08 

5.08 

5.08 

4.82 

4.56 

4.56 

20 

Ca 

6.29 

6.29 

6.31 

5.85 

5.89 

5.89 

22 

Ti 

4.91 

4.91 

4.90 

4.69 

4.78 

4.78 

23 

V 

3.96 

3.96 

4.00 

3.41 

3.48 

3.48 

24 

Cr 

5.64 

5.64 

5.64 

5.09 

5.12 

5.12 

25 

Mn 

5.48 

5.48 

5.39 

4.69 

4.87 

4.87 

26 

Fe 

7.45 

7.45 

7.45 

6.90 

6.93 

6.93 

27 

Co 

4.87 

4.87 

4.92 

4.32 

4.40 

4.40 

28 

Ni 

6.20 

6.20 

6.23 

5.65 

5.71 

5.71 


Table 5; Number of strong lines in the solar and Arcturus spectra. 


Limiting 

depth 

all lines 

Sun 

atomic only 

Arcturus 

all lines atomic only 

lE-4 

10775 

3409 

31192 

4287 

5E-4 

5600 

2426 

19731 

2819 

lE-3 

4291 

2166 

15919 

2559 

lE-2 

1348 

989 

8758 

1590 


Note. — The number of strong lines (atomic and molecular, or atomic only) in the synthetic spectra of the Sun 
and Arcturus in the range between 1.50 and 1.70 gm in dependence of the limiting line depth. 
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damping and HFS differences from line to line. 

Once we had the LOI, we started evaluating them one by one in order of line strength, with the 
strongest line first. For each line we fit the spectral range within 0.8 A around the line center. To 
account for possible wings of nearby strong lines, we calculated the contribution of all lines within 
18 A around each considered LOI. We evaluated two free parameters: radial velocity changes within 
± 0.25 km s“^ (which also accounts for the uncertainty in the central wavelength) and log{gf), 
which is allowed to vary within the following set of rules: 


1. for DRIO and INT, the lines with measured laboratory values were allowed to vary within 
1 sigma of those laboratory measurements, and all other LOI were allowed to vary by -2 to 
+0.75 dex, and 

2. for DR12, the lines with measured laboratory values were allowed to vary within 2 sigma of 
those laboratory measurements, and all other LOI were allowed to vary by -2 to +0.75 dex. 


The reasoning behind the asymmetric limits on the log{gf) variation is motivated by the fact that 
there are bad and missing lines in the base line list. There are a number of strong theoretical lines 
without observed counterparts, and we want these lines to be strongly suppressed, which sets the 
lower limit of —2 dex in log( 5 f/). There are also a number of observed lines without theoretical 
or laboratory counterparts in the base line list, and we do not want to inflate inappropriate very 
weak lines, which sets the upper limit of +0.75 dex in the log{gf). The line list was updated with 
new log{gf) values on the fly during the fitting. Spectral pieces a r ound each LOI were fitted using 
the downhill simplex optimization algorithm by iNelder &: MeadI (|l965l ) . Because the LOIs often 
overlapped, we performed the evaluation iteratively. 


The treatment of HFS changed between the INT line list and the DRI2 line list. Before DR12, 
no consideration was made to force the different components of an adjusted line to scale with each 
other. Thus, if the code required the strongest component of a line with HFS to be adjusted, then 
no change was made for weaker components. For DRI2, we changed HFS component log( 5 f/) values 
by the same amount. 


For DRIO and INT the majority of the LOI were fit only to the Sun, but those lines that were 
not visible in the Sun (i.e., were below the depth threshold) but were visible in Arcturus were fit to 
Arcturus. LOI that were above the detection threshold in the Sun were fit again and superseded 
any corrections that were based on the Arcturus spectrum. Thus, LOI were either fit to Arcturus 
or the Sun. Atomic lines weaker than the cutoff, i.e., not a LOI, in both the Sun and Arcturus were 
left unchanged from the input line list’s values. So for DRIO and INT, the following methodology 
was adopted: 


1. Fit the LOI in the Sun for two iterations for the log( 5 '/) values. 
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2. Fit the LOI in Arcturus for a single iteration for the \og{gf) values. 


3. Fit the LOI in the Sun for two iterations for the log( 5 '/) values. 


4. Fit the strongest LOI in the Sun for damping parameters. 


5. Fit the LOI in the Sun for three hnal iterations of the log( 5 f/) values. 


For DR12 we implemented the use of both the solar and Arcturus spectra to constrain the 
astrophysical log{gf) values. To accomplish this we determined the astrophysical values from 
Arcturus and from the Sun independently instead of serially (as was done for DRIO and INT) 
and then weighted the solutions, with the Sun getting twice the weight as Arcturus because the 
abundances in the solar photosphere are also confirmed by the meteoritic values. Other weighting 
schemes are possible and will be considered for future line lists. 

We noticed that most values of log{gf) settled down quickly after the second iteration over all 
LOI. We also evaluated the damping constant for the strongest lines, applying the same algorithm 
as for the derivation of the log{gf) values but using only the solar spectrum for the comparison. 
Our standard sequence of the line parameter adjustment process was 2 x Elog{gf), 2 x Fidamp, n 
X Elog{gf), where Elog{gf) and Fidamp designate one iteration over the whole LOI list for log(^/) 
and for the damping constants, respectively. 

We justified the number of necessary iterations by running a few long sequences “2 x Elog(gf), 
2 X Fidamp, n x Elog{gfy^ with n over 30. It was noticed that the log(g'/) values for most of the LOIs 
settled down after the “2 x Elog{gf), 2 x Fidamp, 2 x Elog^gf)” sequence. However, the log( 5 f/) 
in a few dozen of the LOIs did not settle down even after 30 iterations. Instead, the parameters 
oscillated around certain values. This usually happened to overlapping lines with similar depth. We 
identified these lines in the last iterations in the “2 x Elog{gf), 2 x Fidamp, 4 x Elog{gf)” sequence 
and replaced these \og{gf) values with their average values. Figure [2] shows a “hair-diagram” for a 
small piece of spectrum with an illustration of log(( 7 /)-stable and unstable lines. Each curve starts 
at its central wavelength at the x-axis. The x-deviation from the vertical designates the difference 
of log( 5 '/) from the original value. The vertical direction shows the iterative sequence progress (2 
X Elog{gf), 2 X Fldamp, and then n x Elog{gf)). 

The version of MOOG we employed to generate the astrophysical line parameters does not have 
a proper treatment for the damping of hydrogen lines. Thus these features were artificially removed 
from the atlases by dividing by a synthetic spectrum containing only hydrogen lines, where these 
broad features were forced to ht. The hydrogen lines were also removed from the input line list 
when importing it into MOOG; the final adopted line list used by APOGEE included the hydrogen 
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1.532 1.533 1.534 1.535 1.536 

Wovelength, mu AND delto(gf) 


Fig. 2.— Example of log{gf) stability during fits to the Sun and Arcturus. The X-axis for this 
diagram is in microns combined with a scaling of offsets in log(gf). The horizontal bar at the top 
right corner corresponds to the length of a 1.0 dex change in log(g'/). Change of the log(g'/) with 
respect to the original value is shown as the deviation from the vertical line. The Y-axis shows 
the progress of the iterations; 2 x Elog{gf), then 2 x Edamp, and then n x Elog{gf). Parameters 
of most of the lines settle down after the 3rd log(g'/)-iteration. The line thickness denote the 
line depth with the thinest lines having depths less than 0.01 and the thickest lines having depths 
greater than 0.2. 
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lines. Molecular lines could be adjusted in a way similar to that adopted for the lines with HFS, 
where lines from the same multiplet or all multiplets could be adjusted by the same amount. We 
tried fitting the extremely weak CO, CN, and OH features in the Sun by eye for DRIO. Given 
the weakness of the lines in the Sun, no attempt to modify the literature molecular log{gf) values 
were made for INT and DR12, so the final adopted line list is a result of an evaluation of only the 
atomic lines, excluding hydrogen. Figure [3] shows an example of the fitted solar spectrum. The 
rms scatter in the difference spectra (lower panel) is reduced from 0.021 with the base line list to 
0.007 with the final DR12 line list. This improvement is seen across the spectrum and is a strong 
indication that this methodology works. The most important test for any line list is to verify that 
it produces reliable abundances results requiring no or little zero point corrections when compared 


with well established abundance trends and values from the literature. We refer to Holtzman et al 


(j2015l l for such a discussion. 


3.2. Astrophysical log(gf) Values 

As described above, the astrophysical log{gf) values were calculated from lab log{gf) values 
(when available), taking into consideration the errors in those measurements. For DRIO and INT 
the astrophysical log( 5 f/) values were allowed to vary within one sigma of the laboratory errors, while 
for DR12 this was expanded to two sigma. Figure 0] shows the difference between the astrophysical 
and laboratory log(g'/) values (top panel) and this same difference divided by the laboratory error 
(bottom panel). For reference, the blue triangles in this figure show how much the log(g'/) values 
would change if we adopted the same criterium used for the literature empirical log(g'/) values, i.e., 
allowed changes of -2 to +0.75 dex. While some of the red squares and black circles in the lower 
panel of Figure 0] do fall at the one and two sigma limits, most are within the error limits and 
centered around zero, which lends confidence that these astrophysical log{gf) values are of good 
accuracy. The exception may be Ti, where the laboratory log{gf) values tend to be larger than 
the astrophysical log{gf) values in both DRIO and DR12. This may be caused by the temperature 
sensitivity of the Ti I lines and thus were driven by the assumed abundance of Ti in Arcturus (the 
cooler of the two stars constraining the astrophysical log{gf) values). The actual errors on the 
laboratory Ti I log{gf) values are very small and therefore the astrophysical \og{gf) values were 
not allowed to vary substantially from the original laboratory log{gf) values. When the laboratory 
constraints were removed for Ti I (the blue triangles in the upper panel of Figure 0|), the resulting 
log(fl'/) values were again not substantially different from the laboratory values, with the single 
exception of the Ti I transition at 16413.029 A. 


3.2.1. Impact of using COD vs. Total Flux 

As mentioned earlier, due to an incorrect entry in the MOOG parameter file, that was rec¬ 
ognized later on, the solar \og{gf) values of the INT and DR12 line lists were calculated in total 
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Fig. 3.— Example of the input and output spectra from the astrophysical adjustment of the line 
parameters. Top panel: observed solar spectrum (dots) and best-fit synthesized spectrum (solid 
curve). The vertical red bars designate wavelengths of lines whose parameters were adjusted. Some 
visible features are not fit because they are molecular features, and for the INT and DR12 line 
lists we made no astrophysical corrections to molecules. Bottom panel: difference between the 
observed and synthetic solar spectra before (lower curve) and after (upper curve) the iterations. 
The upper curve has been shifted by +0.2 along the vertical axis. The weakest features visible in 
the observed spectrum have depth of 0.004 from the continuum and equivalent widths of <0.6mA, 
at the resolution and SNR of APOGEE spectra these features would not be visible. 
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Atomic Number 


Fig. 4.— A comparison of the astrophysical and laboratory \og(gf) values from DRIO and DR12. 
The top panel shows the differences between the log{gf) values for each line and the bottom panel 
shows the difference between the log{gf) values compared to the laboratory error in the log{gf) 
value. In this hgure the black circles represent the atomic lines with lab data from DRIO and have 
been offset by -0.2 in atomic number for better visibility, the red squares represent the lines with 
lab data from DR12 and have been offset by -1-0.2 in atomic number for better visibility, and the 
blue triangles were calculated without any lab constraints using both Arcturus and the Sun with 
limits of -2 and -1-0.75 dex in log( 5 '/). 
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flux rather than as COD, even though the reference solar spectral atlas is COD. Figures [5] and [6] 
show the impact of the difference in the astrophysical \og{gf) values calculated with the proper 
COD option set in MOOG’s parameter file and the “normal” flux option. The top panel shows the 
log(fl'/) differences between COD and total flux for a fit only to the Sun, as was done for INT, and 
the bottom panel shows the differences for a fit using both the Sun and Arcturus, as was done for 
DR12. The Arcturus residuals are best fit with a constant offset of 0.019 dex in both hgures. The 
INT residuals are best fit with a second order polynomial with respect to the excitation potential 
or a constant offset of 0.063 dex as a function of atomic numbers. This shows that by not using the 
appropriate COD flag set in MOOG’s parameter file is more problematic for the INT line list than 
for DR12. The Sun-only fit could introduce a small temperature bias; a simple test using MOOG 
and these two line lists suggests a temperature offset of ~ AOK and a modest metallicity bias of 
~ —0.05 dex, in the sense that the logfqf) values s hould have been larger and thus the derived 
abundances are too large. iGarcfa Perez et al.l (1201311 adopt ed t he INT line li s t and thus was im¬ 
pacted by both of these zero points while ISmith et al.l (j2013l i and iGunha et al.l (j2015l l are impacted 
by only the metallicity bias because their stellar effective temperatures and surface gravities were 
set with an independent methodology. The exact impact these biases have are somewhat dependent 
upon the lines adopted. For example, all abundances based on molecules should not be influenced 
because these lines were not adjusted after DRIO in the astrophysical log{gf) methodology. Lines 
with laboratory constraints may not have been impacted depending on whether the astrophysical 
log(fl'/) derived was at the limits imposed by the laboratory errors. 


In summary, the mean impact of this improper treatment on DR12 is less than 0.02 dex on the 
log(^/) values. In addition, there is no apparent slope with excitation potential or atomic number. 
Since ASPGAP determines stellar parameters using chi squared fits to the entire spectrum, the 
log(fl'/) values for the molecular and hydrogen lines are unchanged by the astrophysical log(g'/) 
corrections, and the errors introduced by using the solar spectrum are, on average, very small, we 
conclude that the impact is well within the errors presented in DR12. 


3.3. Astrophysical Damping Values 

As mentioned in the previous section, after several iterations of the astrophysical log{gf) code 
we fit the line profile for the strongest lines by allowing the lines to be adjusted up to ± 0.2 dex 
from the input damping value. Figure [3 shows the difference between the input and astrophysical 
damping values. The average difference for all damping values is 0.004 dex. Thus, on average, the 
astrophysical code is not making large changes. However, it makes a small systematic change to a 
few elements; for example, among Mg lines the average difference is 0.04 dex (the effect of damping 
is enhanced), and for the Fe lines the difference is —0.04 dex. 
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Fig. 5.— The difference between the astrophysical log(g'/) values calculated with the proper 
center-of-disk (COD) option set in MOOG and the “normal” flux option set in MOOG, shown as 
a function of atomic number. The top panel shows the differences for a fit only to the Sun, as 
was done for INT, and the bottom panel shows the differences for a fit using both the Sun and 
Arcturus, as was done for DR12. The solid line in both panels is the best fit to the data. 
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Fig. 6.— The difference between the astrophysical log(g'/) values calculated with the proper 
center-of-disk (COD) option set in MOOG and the “normal” flux option set in MOOG, shown as a 
function of the lower level excitation potential. The top panel shows the difference for a fit only to 
the Sun, as was done for INT, and the bottom panel shows the differences for a fit using both the 
Sun and Arcturus, as was done for DR12. The solid line in both panels is the best fit to the data. 
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Fig. 7.— The difference between the input literature damping and the astrophysical damping 
values for DR12. For elements with five or more points the median value is shown as a solid square. 
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3.4. Performance of the Astrophysical Parameters 


Over the full range of the line list the use of astrophysical log{gf) values improves the fit in both 
the Sun and Arcturus. For the DR12 line list the rms scatter of the difference between the observed 
spectrum and synthetic fit is 0.010 and 0.027 for the Sun and Arcturus, respectively. Without 
the astrophysical values the rms scatter of the differences were 0.017 and 0.032 for the Sun and 
Arcturus, respectively. The improvement is smaller for Arcturus since the spectrum and residuals 
are dominated by molecular features which are not adjusted in our astrophysical methodology for 
DR12. The evidence of how the line list performs is in its abilit y to derive quant i tative abundances 
which can be compared to literature values for the same stars. iHoltzman et al.l (1201.51 ) shows that 


over the range: T^ff = 3800 — 5300 K, log g = 0.5 — 3.8 dex, [M/H] = —2 —1-0.5 dex, ASPCAP 
delivers Te// to within 91.5 K, log g to within 0.3 dex and absolute abundances of 14 elements of 
0.1 — 0.2 dex accuracy. This scatter is in part due to errors in the line list, but also includes errors 
in carbon isotope ratio mismatch, line spread function modeling errors, NLTE affects, and lack of 
spherical model atmospheres. 


Missing Lines 


The analysis of more than a hundred thousand stars in APOGEE allows us to determine 
where the line list is either missing features or has inaccurate log{gf) values. To construct the list 
of missing lines we start by using the average difference between the best fit synthetic spectra and 
the observed spectra shifted to the rest wavelength for stars of specific effective temperatures and 
gravities. We have chosen to construct three average residual spectra, the first from all giant stars 
with Tgff ~4000 K and log gf ~1.0 cgs dex, the second from all dwarf stars with Tg// ~4800 K 
and log g~4.0 cgs dex, and the last one from all dwarf stars with ~4000 K and log g ~5.0 cgs 
dex. These three residual spectra are shown in Figure [8l The left pane shows the entire spectra 
including the gaps between the APOGEE detectors and the right pane shows a small portion of 
these spectra where a few strong features and several weaker features can be seen in more detail. 
We tabulated all residual features with depth greater than 1% (Table [6j). Eor completeness we also 
in cluded the unidenti f ied fe atures in the Be stars observed by APOGEE as telluric calibrators noted 
bv iGhoinowski et al.l (120141 ). We supplemented this list with the strong lines that were not properly 
m odeled in the analy sis of a set of Eourier Transform Spectrograph (FTS) spectra of standard stars 


m 


Smith et al.l (|2013l h Some of these lines fall between the detectors of the APOGEE instrument. 


Eor the coolest dwarfs we suspect that some of the residual features may be H 2 O and EeH, since 
those lines are known to be detected in M dwarfs in the 77-band. As mentioned in Section [2.1.1.71 
we were unable to find an FeH line list that produced significant improvements for the few M dwarfs 
we tested. Euture work on dwarfs includes MARCS/Turbospectrum modeling on these stars by 
including H 2 O and FeH (jZamora et al.l 1201.51 ). 


In principle, shifting the spectra to their rest wavelengths and averaging in that frame should 
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have eliminated residual features such as badly removed telluric emission and absorption features, 
as well as any interstellar features. However, due to the complicated selection biases in APOGEE 
targeting it is possible that some interstellar or telluric features populate Table [6j One example 
is the feature at 15272.4 A which is the strongest dif fuse interstellar band (DIB) in the H-band 
and is detected in the majority of APOGEE spectra (|Zasowski et al.ll2015l b We hope that future 
versions of the APOGEE line list will include identifications of these types of features. 


5. Summary 


The machine readable tables and data files contain the DRIO, INT an d DR12 line l ists in the 
ASPGAP format, see Tables [71-P Ad ditional copies formatt ed for MOOG (Sneden 1974), Syn spec 


Plea I 2 OI 2 II are 


( Hubenv 2006 : Hubenv k, Lan 2 2011) and Turbospectrum ( Alvarez fc PlezI 199S; 
available in the electronic edition for download. The ASPGAP formatted line lists contain the 
baseline line list, as well as the astrophysical log{gf) and damping values. This line list is formatted 
in a way similar to that employed by Kurucz in his online database and CDROM releases, where 
wavelengths are given vacuum and in nm. Note that over the wide range of this line list, the choice 
of what methodology to use to translate fr om air to vacuum can have a significant impact. We have 
adopted the conversions of ICiddorl (jl996l ) . See Appendix |A] for further discussion on this choice 
and coefficients for conversion. Eor the upper and lower energy levels we use the unit cm~^. 

Some attention has been paid to HES within these line lists but in the very strongest lines of A1 
the lack of the HES components means that the abundances derived are likely to be over-estimated. 
Users of this line list and the abundances of A1 should be wary in the strongest lined and coolest 
stars. 

This line list has been tuned to be applied to red giants over a range that covers E,G,K and 
warm M giants. Tests conducted by the APOGEE team suggest that it also performs adequately 
for F, G, and K dwarfs. These line lists does not contain the FeH lines nor the H 2 O lines needed 
to model cool M giants or M dwarfs. The reader is cautioned when using these line lists outside 
these stellar types or on stars with extreme abundance patterns such as carbon stars where other 
molecular features may play an important role and these line lists have not been tested. In addition, 
the use of the astrophysical log(g'/) values derived for this work should be approached with caution 
w hen using a dif f erent spectral synthesis code or suite of model atmospheres. See Sections 3 and 4 
of IZamora et al.l (|2ni5l ) for more extensive information about the tests that have been conducted 
and possible systematics that could be introduced. 


The ASPGAP machine readable line lists (Tables [7] - [9]) , are 
and ~ 134, 000 rows of data. 


24 MB in size with 35 columns 
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0 


- 0.1 


- 0.2 


-0.3 


-0.4 

Fig. 8.— The average residual spectra (observed spectra - best fit synthetic spectra) for stars near 
Teff ~ 4000 K and log( 5 r) ~ 1.0 cgs dex (top), Te// ~4000 K and log g ~5.0 cgs dex (middle), and 
Teff ~4800 K and log g~4.0 cgs dex (bottom). The left pane shows the entire APOGEE region; 
the breaks in the spectrum are the locations of the gaps between the detectors. The right pane 
contains a small region of these average residual spectra. 
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Table 6: Missing Lines 


A (vac A) 

EW (mA) 

notes 

15174.1 

18 

a 

15177.5 

12 

h 

15178.6 

25 

a 

15180.4 

20 

a 

15182.3 

61 

a, b, h 

15187.3 

19 

a 

15201.1 

16 

a 

15202.4 

24 

a, b, h 

15209.7 

39 

a 

15210.5 

17 

h 

15211.0 

33 

b 

15212.5 

48 

a 

15214.5 

16 

a 

15216.3 

20 

a, b 

15219.4 

59 

b 

15227.6 

10 

h 

15228.7 

12 

h 


Note. — Table [6] is published in its entirety in the electronic edition of the Astrophysical Journal Supplements. A 
portion is shown here for guidance regarding its form and content. 

“Seen in the cool APOGEE giant spectra residual 

*’Seen in the hot APOGEE spectra residual 

“Between APOGEE detectors 

‘'all FTS stars 

“only M FTS stars 

^found only in HD199799 

®see S. Hasselquist et al. 2016 in prep. 

'“Seen in cool APOGEE dwarf spectra residuals 
*Seen in APOGEE Be spectra, see Choinowski et al. ( 2014 1 
^This may be a DIB feature, see lZasowski et al.l (120151 1 
'‘This may be a DIB feature, see Geballe et al. ( 2011 1 
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Table 7; ASPCAP Line List Format for DR12 (20131216) 


Bytes 

Format 

Label 

Description 

1- 9 

F9.4 

lam 

wavelength in vacuum nm 

11- 17 

F7.3 

orggf 

original log{gf) value 

19- 25 

F7.3 

newgf 

improved literature or laboratory log{gf) 

27- 30 

F4.2 

enewgf 

error on \og{gf), when available 

32- 34 

A3 

snewgf 

source for improved literature or laboratory log(gf) 

35- 41 

F7.3 

astgf 

astrophysical log{gf) 

43- 45 

A3 

sastgf 

source for astrophysical log{gf) 

47- 54 

F8.2 

specid 

species id 

55- 66 

F12.3 

EPl 

lower Energy Level in cm“^ 

67- 71 

F5.1 

J1 

J value for EPl 

72- 82 

All 

EPlid 

EPl level identification 

83- 94 

F12.3 

EP2 

upper Energy Level in cm“^ 

95- 99 

F5.1 

J2 

J value for EP2 

100-110 

All 

EP2id 

EP2 level identification 

111-116 

F6.2 

Rad 

Damping Rad 

117-122 

F6.2 

Sta 

Damping Stark 

123-128 

F6.2 

vdW 

Damping vdW 

130-131 

12 

unite 

NLTE level number upper 

132-133 

12 

Inlte 

NLTE level number lower 

134-136 

13 

isol 

isotope number 

137-142 

F6.3 

hyp 

hyperfine component log fractional strength 

143-145 

13 

iso2 

isotope number 

146-151 

F6.3 

isof 

log isotopic abundance fraction 

152-156 

15 

hEl 

hyperfine shift for first level to be added to El in mK 

157-161 

15 

hE2 

hyperfine shift for first level to be added to E2 in mK 

162 

A1 

EO 

hyperfine F symbol 

163 

11 

El 

hyperfine F for the first level 

164 

A1 

notel 

note on character of hyperfine data for first level: z none, ? guessed 

165 

A1 

S 

the symbol for legibility 

166 

11 

E2 

hyperfine F’ for the second level 

167 

A1 

note2 

note on character of hyperfine data for second level: z none, ? guessed 

168-172 

15 

gl 

lande g for first level times 1000 

173-177 

15 

g2 

lande g for second level times 1000 

178-180 

A3 

vdWorg 

source for the original vdW damping 

181-186 

F6.2 

vdWast 

astrophysical vdW damping 


Note. — 

In addition the same DR data is available in MOOG, Tnrbospectrum and ASSeT formats in the included .tar.gz 
package. 

This table is available in its entirety in FITS format. 

This table is available in its entirety in machine-readable form. 
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Table 8: ASPCAP Line List Format for INT (20120216) 


Bytes 

Format 

Label 

Description 


Same Format as Table [7] 


Note. — 

In addition the same DR data is available in MOOG, Tnrbospectrum and ASSeT formats in the included .tar.gz 
package. 

This table is available in its entirety in FITS format. 

This table is available in its entirety in machine-readable form. 









- 31 ~ 


Table 9: ASPCAP Line List Format for DRIO (20110510) 


Bytes 

Format 

Label 

Description 


Same Format as Table [7] 


Note. — 

In addition the same DR data is available in MOOG, Tnrbospectrum and ASSeT formats in the included .tar.gz 
package. 

This table is available in its entirety in FITS format. 

This table is available in its entirety in machine-readable form. 
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For a detailed discussion on how these line lists pe rformed we point th e reader to lMeszaros et al 


( 201311 and Ilya, mb et al. ( 2015) for DRI P performance, Smith et al.l ( 20131 ) for discussi ons of the INT 
l ine lis t and iMeszaros et alJ (|2015l ) and iHoltzman et alJ (j2015l ) for DR12 performance. iMeszaros et al 
(|2013l ) compared literature metallicity, effective temperature, and surface gravity to the values com¬ 
puted by ASPCAP for DRIO. They found some systematic offsets in gravity and metallicity. The 
source of the se offsets ma. v not b e the line list but rather the methodology within the DRIO version 
of ASPCAP. iLamb et al.l ()2015l ) conducted an independent manual analysis of several metal-poor 
globular clusters stars with optical spectra and APOGEE DRIO spectra and line list. They found 
good agreement for Ee, Mg, and C a lines, but had a concern about the single detectable Ti I 
line (15339 A vac) being a blend. iLamb et al.l ( 20151 ) adopted an average between the optical 
and ff-band result s for Fe, Mg, and Ca and adopted the i^-band results for Al, Si, O, C, and N. 
Smith et al.l (j2013l ) conducted an independent manual analysis of several very high resolution, very 
high S/N FTS spectra using the INT line list. In a comparison of their results to the literature, 
they found agreement within ~0.1 dex for a ll abundances they d erived: ^^C, ^^C, N, O, Mg, Al, Si, 
K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, and Cu. ICunha et al.l (j2015l ) analyzed manually, also using the 
INT line list, a sample of red giants and clump stars in the very metal rich open cluster NGC6791 
and found overall good agreement with optical results from the literature for the abundances of 
oxygen, sodium and iron. 


Meszaros et al.l (j2015l ) derived abundances in an independent manual analysis of APOGEE 


spectra of 10 globular clusters using the DR12 line list. They derived abundances for nine ele¬ 
ments and found good agreement for most elements, once zero point differences in adopted solar 
abundances are considered. Exceptions included Ti and Ca, which exhibited a large scatter in the 
APOGEE results, and Al, which are based on ver y weak spectral lines in the literature optical 
analyses and on very strong features in the 77-band. iHoltzman et al.l (1201.51 ) analyzed the ASPCAP 
DR12 results for self consistency within clusters and also did comparisons with abundance results 
in the literature. They found that DR12 results had an internal abundance consistency at the 
level of 0.05 - 0.09 dex and 0.1 - 0.2 dex agreement with literature values. This study, however, 
pointed out a number of elements which exhibited unexpected trends with respect to metallicity. In 
particular, for the entire APOGEE sample, the mean abundances of S, Si, and Ca at roughly solar 
metallicities were above the solar value, and Ti does not exhibit the expected rise at decreasing 
metallicity seen in the literature. These types of analyses and, in particular, comparisons with 
cluster abundance results can act as a guide for future improvements to the line list. 


A. Air-Vacuum Conversion 

Since the APOGEE spectrograph is operating in vacuum, it was natural for APOGEE to 
adopt a vacuum wavelength scale. Further support for embracing a vacuum scale is motivated by 
the fact that, historically, the definition of standard air includes tight constrains on temperature 
(15 C), pressure (1 atmosphere), and humidity (dry), but not so clearly the C02 concentration, 
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which is time and spatially dependent, and the standard temperature scale has also seen multiple 
changes in the past century, and these drag the air-to-vacuum corrections. In addition, the reference 
wavelengths for the Th and Ar lines usually employed for calib ration come from vacuum-housed 


spectrographs, and are therefore more naturally used in vacuum ( 

Morlra 

1973; 

Palmer &: Eneleman 

1983: 

Hinkle et al. 

2001; 

Lovis &: Pepe 

2007: 

Kerber et al. 

200f 

1). This appendix describes the 


differences on the vacuum-to-air conversions and motivates APOGEE adopting a specific one for 
dealing with such corrections. 


A.l. Available formulae 


The lAU standard for the vacuum to standard ai r corrections (see resolution No. C15, Com- 
mi ssion 44, XX I General Assembly in 1991) refers to[Oosterhofl (jl957l ). which adopts the results 
by Edlen (1953). 


Ao-A 


= n — 1 = a + 


61 


+ 


62 


(Al) 


A " ' cl — 1 /Aq c2 —1/Aq 

where Aq is given in /rm, and the cons tants are give n in the hrst row of Table [TOl Later work widely 
adopted in the phy sics literature by lEdlenl (119661 ) rederived the constants from optical and near 
UV data, and later iPeck fc Reeder! (119721 . and references therein) added additional measurements 
extending into the IR, up to 1.7 /rm, which showed a systematic deviation from Edl&s equation at 
the level of several 10“®, or a few ms“^. 


Table m compares t he parameters proposed by lEdlenl ([ 


(11972), and 


Giddod (jl996l ) , the latter reference largely based on 


(1953) 

Edlen ( 

196 

6), 

Peck Kr. Reeder 

a Peck & E eeder (' 

. 972 IL but updated 


to account for the changes taken place in the international temperature scale since their work, and 
adjusting the results for the CO 2 concentration. Fig. [9] shows a maximum disc repa ncies at the level 
of 3 X 10“®, or 10 ms“^ at 1.6 pm. The equations proposed bv lEdlenl (jl966l j and iPeck &: Reeder 


(jl972 ) differ only by about one fourth of that. About half of the difference between Peck & Reeder 
and Ciddor (2 x 10“® or 6.5 ms“^ at 1.6 pm) is related to the temperature scale. The paper by 
Peck &: Reeder is not specific regarding the scale used, and Ciddor assumes they used the IPTS-48 


Table 10: Parameters for Eq. lAll 


Reference 

a 

bl 

b2 

cl 

c2 

Edl& 119531 

6.4328x10"® 

2.94981x10-2 

2.5540x10"^ 

146.0 

41.0 

Edlen 119661 

6.4328x10"® 

2.406030x10-2 

1.5997x10"^ 

130.0 

38.9 

Peck & Reeder 119721 

0.0 

5.791817x10-2 

1.67909x10"® 

238.0185 

57.362 

Ciddor 11996) 

0.0 

5.792105x10-2 

1.67917x10"® 

238.0185 

57.362 
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standarci^, which is warmer than the one currently in use (ITS-90) by 9.2 mK at 15 C. If Peck & 
Reeder used the IPTS-68 standard instead, the difference would be reduced to 5.6 mK. Changes of 
9.2 mK and 5.6 mK amount to an increase in (n — 1) xc of about 2.6 and 1.6 m“^, respectively, 
at 1.6 fiia. The rest of the correction between Peck &: Reeder and Ciddors review (3.9 ms“^ at 
1.6 fj,m) is related to the CO 2 concentration in standard air. Because of secular variations in the 
typical laboratory air, Birch & Downs proposed to use 450 ppm in the definition of standard air, 
while Ciddor estimates that a value closer to 300 ppm is adeq uate for the me a surein ents used by 
Peck fc: Reeder. In summary, the variation between between iPeck &: Reeder! (jl972l l and [Ciddor 
(| 19961 1 is connected to a change between the actual conditions of standard air now and at the time 
of the measurements. 


A.2. Conclusions 


In view of the preceding discussion, we underline the advantages of using vacuum wavelengths 
for the APOGEE spectra. In some cases, corrections bet ween air and v acuum will be needed, and 
for those it is proposed that the formulation proposed bv ICiddon (jl996l ) be used. This corresponds 
to Eq. lAll with the Ciddor constants given in Table m This is valid for the wavelength range 
between 0.23 and 1.7 fiui, and the estimated accuracy in the predicted refraction index of standard 
air is about 10“® or roughly 3 ms~^ at 1.6 /rm. Note that it is straightforward to go from vacuum 
to standard air wavelengths with Eq. \AT\ but the inverse process requires iteration since n is given 
as a function of vacuum wavelength. 
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^ Intern ational Practical Temper ature Scale; tes 
(100 — t 4 s) i Fofonoff fc Brvden 1975 1 


= 1.00024 X tgo ( Saunders 1990h : tes = tis — 4.4 x 10 ®t48 
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Wavelength {fu.m) 


Fig. 9.— Difference between the refractive index of standard air (n) and unity in the ff-band for 
four different sources considered. 
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